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ABSTRACT 


The  interaction  of  internal  waves  and  surface  waves  in 
water  is  explored  in  the  regions  where  the  effects  of  the 
interaction  are  small.  Explicit  solutions  for  an  arbitrary 
internal  wave  are  given  for  various  kinds  of  boundary  con¬ 
ditions,  with  a  detailed  discussion  of  the  conditions  under 
which  they  are  valid.  The  connection  of  the  solutions  with 
the  conservation  equations  is  made  explicit,  and  the  modifi 
cations  produced  by  viscous  damping  are  outlined. 
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I.  INTRODUCTION 


This  paper  is  to  be  viewed  as  a  sequel  to  an  earlier  work 
on  the  internal  wave-surface  wave  interaction  in  water,  which 
appeared  in  IDA/ JASON  Study  S-334.*  We  wish  to  enlarge  on  the 
previous  analysis  in  the  following  ways: 

1.  Clarify  the  assumptions  involved  in  the  calculation. 

2.  Correct  the  result  in  (I)  by  including  a  term  which 
was  unjustifiably  neglected. 

3.  extend  the  solution  to  ranges  of  the  parameters  not 
treated  in  (I). 

4.  Establish  the  connection  with  the  approach  to  the 
problem  based  on  the  conservation  equations,  and 

5.  Include  the  effect  of  viscous  damping  of  the  surface 
waves . 

Let  us  begin  with  a  brief  recapitulation  of  the  problem 

which  we  wish  to  discuss.  We  assume  the  existence  originally 

of  a  surface  wave,  that  is,  of  a  displacement  of  the  surface 

of  an  infinite  ocean,  given  by  h  (x,t)  =  AQ  exp  i(kQx  -  a>Qt) 

with  u>  =  /g'F  . **  An  internal  wave  with  velocity  U- -which  is 
o  °  o 


The  complete  reference  is  IDA/JASON  Study  S-334,  Generation 
and  Airborne  Detection  of  Internal  Waves  from  an  Object  Moving 
Through  a  Stratified  Ocean,  JASOhl  1968  summer  study,  Vol.  1 1 , 
April  1969,  p.  69.  We  shall  refer  to  this  as  Cl)  in  what 
follows . 

We  shall,  for  simplicity,  assume  there  is  only  one  horizontal 
dimension  (denoted  xj .  The  depth  variable  is  called  z,  and 
is  measured  positive  up.  We  are  also  going  to  confine  our¬ 
selves  to  gravity  waves- -hence  (since  the  water  is  assumed 
to  be  infinitely  deep)  we  take  u)Q  =  /gkQ .  This  limits  us  to 
values  of  1;  <  3  cm’1.  ° 
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assumed  to  be  &  given,  but  arbitrary,  function  of  space  and 
time- -is  turned  on,  and  this  produces,  for  later  times,  a 
change  6h(x,t)  in  the  surface  displacement.  The  problem  is  to 
calculate  6h. 


We  ..hr7'.  I  carry  out  the  calculation  under  the  following 
as  sump  ■.  i  O’u* 

1.  '^mall  amplitude  surface  waves: 

kA  «  1  (1.1) 

2.  Small  internal  wave  velocities: 

U/c  «  1  (1.2) 

o 

where  c^  =  1_  /g/k  is  the  group  velocity  of  the 
surface  wave. 

3.  Long  wavelength  and  long  period  internal  waves: 

X  <<  k  and  9  «  to  (1.3) 

where 

„  „  1  3U  and  „  n  1  3U 

K  '  u  n  n  u  at 


are  the  wave  number  anu  frequency  of  the  internal  wave. 

4.  The  internal  wave  effects  on  the  surface  wave  are 
smal i : 

«  1,  Sk/k  <<  1,  jp  «  1  (1.4) 

Thus  a  perturbation  treatment  of  the  effect  is 
appropriate . 

5.  In  (I),  we  also  made  the  assumption  that  the  time  of 
interaction  between  the  surface  wave  and  internal 
wave  was  short: 

(K/k)  (uit)  «  1  (1.5) 

If  viscosity  damps  the  surface  wave  in  a  time  1/a, 
this  condition  is  satisfied  if  (K/k)(ai/o)  <<  1. 
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Assumption  (5)  was  made  for  convenience  in  (I),  but  was 
not  essential  to  the  t; ^atnent .  We  shall  therefore,  here, 
also  obtain  the  solution  without  imposing  this  constraint.  We 
shall,  however,  s.ill  find  it  useful  in  practice  to  keep  a 
weaker  limit  on  the  "ime,  namely  (K/k)4(u>t)2  <<  1;  to  relax 
this  requirement  too  is  possible,  but  prevents  us  from  exhibit 
ing  the  solution  in  a  very  explicit  form. 

A  typical  set  of  numbers  with  which  we  may  be  concerned 
is  as  follows: 

k  =  1  cm" 1 

A  =  0.1  cm 

U  =  1  cm/sec 

K  =  10’ 4  cm’ 1 

C  =  15  cm/sec 

For  k  =  1  cm’1,  we  have  a  =  3x10“ /sec,  c  =  15  cm/sec  and 

o 

u  =  50  sec’ 1 . 

With  these  parameters  all  of  our  conditions  would  appear 
to  be  met:  we  have 

kA  =  0.1 

U/c  =  0.07 
8 

K/k  =  10’4 

n/w  =  5x10’ s 

(K/k)  •  a i/a  =  0.1 

The  question  of  how  well  this  idealized  problem  imitates 
the  real  ocean  is  of  course  a  serious  one.  Probably  the  most 
important  physical  effect  that  has  been  ignored  entirely  is 
wind.  We  have  pretended  that  the  wind  enters  only  indirectly, 
in  that  it  is  responsible  for  the  generation  of  the  surface 
waves,  but  that  once  generated,  the  surface  wave  is  no  longer 
seriously  influenced  by  the  wind  and  is  merely  eroded  by 
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viscosity  in  some  finite  distance- -or  time- -depending  on  its 
wavelength.  This  is  evidently  a  great  idealization. 

The  next  most  serious  difficulty  would  appear  to  be  in 
the  validity  of  the  kA  <<  1  assumption.  Nonlinear  effects  in 
the  surface  wave  amplitude  may  well  become  important  as  SA/A 
increases  from  zero,  and  invalidate  the  treatment  given  here 
even  when  6A/A  is  still  quite  small. 
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II.  BASIC  EQUATIONS 

To  obtain  our  basic  equations  we  first  write  the  equations 
for  a  free  surface: 

9E  8  0  BIT  _  5  4> 

Jt  +  Hx  9x  =  9  z 


i£  ♦  1 

9t  7 


•V  + 

c  J  \9  zj 


=  -gH  , 


both  evaluated  at  z  =  E.  The  notation  is  obvious:  E  is  the 
surface  displacement  and?  is  the  velocity  potential.  We  have, 
in  addition 

V2?  =  0 

and  the  boundary  condition  that  4>  -*■  $ ,  where  7<J>  =  U,  as 
z  -*•  -  <I>  is  assumed  to  be  a  given  function,  and  represents 

the  velocity  potential  of  the  internal  wave. 

The  next  step  is  to  expand  the  surface  equations  to  second 
order  in  E.  Then  we  write 

E  =  h  +  H 

<p  =  <p  +  $ 

and  keep  all  first-order  term;  in  either  h  or  <t> .  We  choose  H 
and  *  (which  are  of  course  to  be  identified  with  the  internal 
wave)  to  satisfy  the  free  surface  equations  by  themselves  (to 
fust  order  in  H--and  we  furthermore  assume  that  the  z-depend- 
ence  in  the  internal  wave  is  weak.)  We  then  obtain  equations 
for  h  and  <P ,  viz: 


9h  A  „  9h  .  .  9U  _  9<f>  .  1  9£  311 
9t  +  U  97  h  9x  "  9z  g  9x  9t 
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(2.1) 


i 


3(J)  9<J> 

Jt  +  u  lx  =  -gh 

evaluated  at  z  =  0,  together  with 


(2.2) 


V24>  =  0 

and  the  boundary  condition  <j>  ->  0  as  z  +  •  »,  The  symbol  U 
stands  for  the  x  component  of  the  internal  wave  velocity  ft. 

We  have  thus  far  made  use  of  assumption  1.1,  the  small 
amplitude  assumption. 

Equations  2.1  and  2.2  are  our  starting  point;  comparing 
with  (I)  we  note  the  presence  in  Eq.  2.1  of  the  term  h  9U/9x. 

In  (.1)  this  was  neglected  in  accord  with  assumption  1.3  with 
the  argument  that  it  was  small  compared  to  U  9h/9x.  In  mag¬ 
nitude,  this  is  true;  but  the  two  terms  are  of  different  phase, 
and  hence  it  is  not  legitimate  to  ignore  h  3U/9x.  This  fact 
will  become  more  transparent  later. 

The  next  convenient  thing  to  do  is  to  eliminate  <p  between 
Eqs .  2.1  and  2.2  to  obtain  an  equation  for  h  alone,  which  is, 
after  all,  the  quantity  of  primary  interest.  This  is  made 
feasible  through  V2<j)  =  0  and  the  fact  that  <J>  0  as  z  -  °°, 

which  permits  us  to  write 

,  ,  .  r  dk  *  ^  ikx  I  k  I  z 

<J>(x,z,t)  -  f  -yr-  4>(fc,t)  e  e 


Hence  (we  take  k  positive)  we  have 

9<J>/9z  =  -i  34>/9x  (2.3) 

We  can  now  get  rid  of  <j>.  Equations  2.1  and  2.2,  together  with 
assumptions  1.2  and  1.3,  yield 


•  3  \  v.  on  3Zh  7  9U  3h  7  3U  9h 

Tx  *  ^  ~  -  ‘  L  TT  TTT  L  rT 


TficlTt  '  It  7x  lx  It 


(2.4) 


We  have  neglected  all  quadratic  terms  in  U,  and  all  second  and 
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higher  derivatives  o£  U,  in  obtaining  this  result.  The 
boundary  condition  which  we  associate  with  Eq.  2.4  is  that  for 
t  >  0,  h  =  hQ(x,t),  some  specified  surface  wave  present  before 
the  internal  wave  was  turned  on.  Normally,  we  shall  take  hQ 
to  be  simply 

h  =  A  e^V'VJ  (2.5) 

0  0 


where 


pF-j\V^Hijfc’J*V(*7^)  -*»<*  WJi 


TTfe^TS^.T?  aj"?;  ^»* 


III.  PERTURBATION  SOLUTION 

We  assume  that  the  effect  of  the  internal  wave  on  the 
surface  wave  which  was  present  initially  will  be  small  (As¬ 
sumption  1.4).  We  may  then  write 


h  =  h  +  6h,  h  >>  6h 
o  o 


(3.1) 


where,  to  first  order  in  the  perturbation, 


32_ 
^t2 


3  h„  an  3hrt 
ig  ~i  6h  =  -2U  -  2  JiJT 


an  3^ 

9  3U  O 

2  3x  "5T  • 


for  the  simple  choice  (2.5),  this  simplifies  to 

’’  -  ie  y  {h  ■  l"2ko“o  U  *  Zi“o  §  •  2iko  It 


.at 


Z(x,t)  hQ  (x , t) 


(5.2) 


This  will  be  recognized  as  the  same  equation  as  that  used  in 
(I)  except  that  the  coefficient  of  3U/3x  in  (I)  was  iuQ  rather 
than  2ici)Q.  The  factor  of  2  reflects  the  existence  of  the 
h  3U/3x  term  in  Eq.  2.1. 

It  may  be  objected  that  the  derivation  of  Eq.  3.2  is  in¬ 
valid,  since  it  relied  on  the  assumption  that  hQ  >>  6h,  and 
for  certain  values  of  x  and  t  (for  which  hQ  vanishes)  this  in¬ 
equality  is  false.  However,  these  values  constitute  a  very 
small  range  of  x  and  t,  and  do  not,  in  fact,  affect  the  cor¬ 
rectness  of  Eq.  3.2.  This  will  be  made  explicit  in  Section  IV, 
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where  we  will  obtain  the  same  result  without  use  of  the  as¬ 
sumption  that  hQ  >>  6h. 

In  any  event,  our  problem  is  now  a  straightforward  one, 
namely  to  solve  Eq.  3.2  with  the  boundary  condition  that 
6h  =  0  for  t  <  0.  This  can  be  done  directly  by  the  method  used 
in  (I).  Let  A  be  the  retarded  Green’s  function,  so  that 

fe  -  i“  fe)4  -  «00  «(t)  (3.3) 

with  A  =  0  for  t  <  0. 

Then 


6h(x,t)  =  //  dx'dt'  A(x-x’,t-t')  Z(x’,t')  h0(x’,t') 


and  using  the  (obvious)  solution  to  Eq.  5.3  together  with 
Eq.  2.5,  we  find 


6h (x ,  t) 
hQ(x,t) 


-  /  Jf  I 


dk 

2tt 


Z(k.t') 


„ikx 

e 


dto  e 


-  iu)  Ct- 1 ' ) 


Tit  gk-  2wu)o  -to5 


(5.4) 


where 


Z (k, t)  =  /  dx  e* lkx  Z (x , t)  . 

"  00 

The  integral  over  eta  in  Eq.  3.4  is  to  be  taken  above  the  poles 
of  the  denominator  in  the  complex  w  plane  because  A  is  the  re¬ 
tarded  Green's  function. 

Up  to  this  point  we  have  invoked  assumptions  1.1,  1.2, 

1.3  and  1.4  but  have  made  no  use  of  1.5.  Equation  3.4  con¬ 
stitutes  a  complete  solution  to  the  problem;  however  for  a 
given  (arbitrary)  U  some  numerical  integrations  will  be  re¬ 
quired  to  obtain  a  solution.  It  is  therefore  of  some  interest 
to  try  to  simplify  3.4. 
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To  do  this  we  note  that  from  assumption  1.3,  U,  and  hence 
Z,  may  be  expected  to  have  only  small  frequency  components 
present.  Hence  we  may  expect  ua  to  be  small,  so  we  may  expand 
the  denominator  in  Eq.  3.4: 


<5  O  >  / 

gk-2 uiu^-ur  gk-2u)too  (gk- 2!i)mo) 2 

The  first  term  in  this  expansion,  after  evaluation  of  the  con¬ 
tour  integral  over  du),  yields  directly 

ir  -  dr  Jl  dt'  2(*  -  M'-1')*  »')  cs-*) 

O  O  -OO  6 

Lot  us  assess  the  error  involved  in  stopping  at  this  point. 

To  do  this  it  is  important  to  note  that,  from  Eq.  3.2,  Re  Z  is 
bigger  than  Im  Z  by  a  factor  'v  k/K.  Since  the  next  term  in 
the  expansion  is  out  of  phase  by  90°  with  this  term,  the  first 
term  alone  will  be  an  excellent  approximation  to  Im  6 h / h q  but 
a  less  good  or.e  for  Re  oh/hQ.  In  fact  the  errors  are  evidently 
of  order  0Qt)  (K/ko) 3  and  (u)Qt)  (K/k^)  for  the  imaginary  and 
real  parts,  respectively. 

The  first  term  is  thus  a  good  approximation  for  the  real 
part  of  <$h/hQ;  that  is,  for  &A/AQ,  if  the  "short  time  assump¬ 
tion,"  1.5,  is  valid.  It  is  a  good  approximation  for  the 
imaginary  part,  that  is  for  6w  and  6k,  in  any  event. 

In  the  case  of  short  wavelength  internal  waves,  or  long 
wavelength  surface  waves,  where  assumption  1.5  may  not  be 
satisfied,  we  may  keep  the  next  term  in  the  expansion  3.5. 

This  leads  to  the  result 


10 


■yVTRIT'^'V"'’ 
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r  3  ~ vwtjk '-jy,v  ? 


IT  =  7S~  /  z(x  ’  cg  (t-t'),  t'J  dt' 
o  o  ® 


— •/  (t-fj 


4  03  ‘  -«0 

o 


at,2  z(x,>t,)  •x'=x-cg(t-f)  dt' 


(3.7) 


We  may  again  estimate  the  error  involved  in  stopping  here, 
and  because  of  the  alteration  in  phase  by  90°  of  successive 
terms  in  the  expansion,  it  is  clear  that  this  is  now  an  ex¬ 
cellent  approximation  for  both  real  and  imaginary  parts. 

The  successive  terms  are  of  the  following  magnitudes: 


Imaginary  part:  £  ♦  (^0  *  (^)‘  *  "  •] 
Real  part:  [fe)  *  (^)  *  - '  ] 


Thus  the  neglected  terms  are  of  order  (u)Qt)2  (ft/u)o)4  in  both 

cases . 


Finally,  we  may 
<5A/Aq  and  6k/kQ,  the 
using  the  definition 


extract  from  Eq.  .7  expressions 
quantities  of  real  interest.  We 
of  Z, 


for 
f  ind , 


and 


-  /‘sr  u( 

-  OO 


X 


Cg  £t-t’),  t')  dt' 


t 

J 


-  OO 


(t-t») 


a2 

3x2 


c 

g 


U(x 


(t-t’),  t’)  dt’ 

-  cg  (t-t’),  t')  dt’ 


(3.8) 


(3.9) 
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The  last  term  in  3.9  may  be  dropped  if  the  short  time  assump¬ 
tion  is  valid;  we  have  neglected  the  corresponding  term  in 
<$k/kQ  in  any  event, 

Comparing  3.8  and  3.9  with  the  results  given  in  (I),  we 
note  that  the  expression  for  6k  is  the  same,  but  that  for 
6A  differs  by  factors  of  2  (apart  from  the  32U/3x2  term);  these 
changes  are  caused  by  the  h  3U/3x  term  in  the  original  equation 
with  which  we  began. 
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IV. 


CONSERVATION  EQUATIONS 


We  now  return  to  our  starting  point,  Eq.  2.4.  Instead  of 
making  a  perturbation  expansion  of  this  immediately,  as  we  did 
in  Section  III,  let  us  use  it  to  derive  equations  for  the  ampli¬ 
tude,  wave  number  and  frequency  directly.  This  is  easily 
done  with  the  following  definitions: 


h  - 

AeiX 

(4.1) 

k  = 

3X/3x 

(4.2) 

w  = 

- 3X/3 t 

(4.3) 

The  definitions  4.2  and  4.3  lead  immediately  to  the  so-called 
"conservation  of  waves": 


3k 

Tt 


+ 


3(jJ 

Tx  3 


0 


(4.4) 


We  substitute  4,1  to  4.3  into  2.4,  and  separate  2.4  into 
real  and  imaginary  parts,  and  find  after  a  small  amount  of 
algebraic  juggling  that 

w2  .  gk  -  1  ilA  *  2koi  U 
A  3 1 2 


+ 


2 


3U  1  3A 

Jt'KTx 


3U  1  3A 
Tx  S  Tt 


* 2U  (x  5sfe) 


(4.5) 


and 


Ml  A  3A  .  it  1  9A  .  u 

o)A"5T  £  ITx  ’  oj  ITx 

+  k  au  3U  (4.6) 

a)  Tt  ’  "?x 

The  first  of  these  is  the  real  part  of  2.4,  the  second  the 
imaginary  part. 

These  two  equations  are  recognizable  as  the  dispersion 
relation  and  the  conservation  of  energy,  respectively.  (One 
should  multiply  the  second  through  by  A2  to  make  this  obvious.) 

Let  us  now  introduce  the  perturbation  expansion  in  U. 

Thus  we  write  A  =  A  +  <5A,  k  =  k  +  6k,  u>  =  uj  +  6uj  and  assume 

o’  o  o 

6A  <<  Aq,  etc.  The  zero-order  equations  obviously  tell  us 

that  A^  =  const,  =  /gk  and  k  =  const.  The  first-order 
o  o  °  o  o 

equations  are 

6cj  =  k  U  +  c  6k  (4.7) 

o  g 

(ft  *  cg  f*)sl"  -ko  K  <4'8> 

(l_  +  c  3  5k  =  3  au  +  1  3U  (4.9) 

\  Tt  g  3x  )  IT  4k  cix  ?  5x  4c  Tt 

\  '  oo  g 

where  we  define  c^  =  g/2oJo  =  (1/2) /g/kQ  . 

We  have  used  4.4  to  eliminate  6w  in  deriving  these,  and 
have  also  neglected  higher  derivatives  of  6A. 

First  of  all,  these  equations  are  immediately  seen  to 
coincide  with  Eq.  3.2,  thereby  confirming  the  use  of  hQ  >>  6h, 
in  spite  of  the  small  range  of  x  and  t  over  which  this  may 
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* &WVSZ&* -fc$#WU  •*,  ',-^’*g*i5T? 


*;  5?  “*  r'-VC^^^r-JTVJ ' 


i 


not  be  true.  Second,  these  coincide  with  the  usual  results 
obtained  from  the  conservation  equations.* 

Finally,  the  solutions  may  be  written  down  by  inspection, 
(incorporating  the  boundary  conditions  that  6A  and  6k  vanish 
when  t  <  0)  and  are  seen  to  coincide  exactly  with  Eqs.  3.8  and 
3.9.  The  fact  that  the  solutions  to  4.8  and  4.9  agree  exactly 
with  3.8  and  3.9,  which  are  only  the  first  two  terms  in  the 
expansion  of  the  solution  Eq.  3.4  to  3.2,  stems  from  the 
neglect  of  higher  derivatives  of  6A  in  deriving  4.8  and  4.9. 
This  observation  should  help  to  clarify  what  is  being  dropped 
when  higher  terms  in  the  expansion  3.5  are  neglected. 


0.  M.  Phillips,  Dynamics  of  the  Upper  Ocean,  Cambridge,  1966, 
p.  50.  Noce  that  we  have  an  additional  term  in  9U/3t  not 
given  by  Phillips.  Presumably  he  omits  this  because  of  the 
time  averaging  he  uses  to  obtain  the  conservation  equations. 
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V.  PROPERTIES  OF  THE  SOLUTIONS 

Let  us  discuss  some  special  cases.  First,  suppose  the 
internal  wave  has  a  constant  phase  speed  C,  and  is  turned 
on  at  some  time  t  ,  so  that  U(x,t)  =  U(x  -  Ct)  9(t  -  tQ) . 
Then 

6k  1  fx'Ct  3U(y)  , 

o  g  C  x-cgt+(cg-C)t0  y 

U(x  -  Ct)  -  U(x  -  c  t  +  (c  -  C)  t  ) 


Furthermore, 


k  U  +  c  6k 
°  g 


Hence 


-  C  6k  =  Ufx  -  c  t  +  (c  -  C)  tj 

C-C  ^  2  ^  g  0J 


Thus  if  t  +  *  »,  6u  -  C  6k  -*■  0.  Hence  infinitely  long  after 
the  turn-on  time,  the  quantity  w  -  Ck  is  a  constant.  This,  of 
course,  is  just  the  frequency  in  a  system  moving  with  the  in¬ 
ternal  wave  phase  speed,  and  the  statement  that  this  is  con¬ 
stant  is  part  of  the  "steady- state  solution"  to  the  internal 
wave-surface  wave  interaction  problem. 
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*  t 


The  fact  that  at  finite  times  after  turn  on,  the  steady- 
state  solution  is  not  reached  is  because  the  presence  of 
9(t  -  t  )  in  the  form  chosen  for  1J  violates  the  assumption  of 
the  steady- state  theory,  that  U  is  only  a  function  of  x  -  Ct. 
It  would  thus  not  seem  to  be  possible  to  reconcile  the  steady- 
state  case  with  an  initial  value  problem. 

Next  suppose  we  go  to  resonance,  that  is,  we  take  C  =  c  . 

S 

Then  (choosing  t  =  0  for  convenience) 


and  similarly, 
6A  U 


9U  t2cg  a2u 

3x  "  8  3x2 


(5.2) 


where  we've  kept  the  last  term  in  6a/Aq  in  case  the  "short  time" 
approximation  is  invalid. 

These  results  constitute  the  expression  of  resonant  en¬ 
hancement-  -  they  depend  only  on  the  form  U  =  U(x  -  Ct)  9(t)  and 
the  choice  C  =  c  .  The  linear  (or  quadratic)  growth  with  t 
obviously  means  that  the  expressions  break  down  eventually, 
in  that  the  requirements  6k/kQ  and  6A/AQ  <<  1  fail. 


I 

I 

I 

I 


•-<* 
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VI.  OTHER  BOUNDARY  CONDITIONS 


i 


\ 

1 


Experiments  on  internal  wave- surface  wave  interactions 
tend  to  replace  the  initial  value  problems  we  have  dealt  with 
up  to  now  with  a  boundary  condition  specifying  the  surface 
wave  at  a  point  in  space,  the  point  where  the  surface  wave  is 
generated.  If  it  is  generated  with  constant  frequency,  these 
conditions  would  say  that  at  x  =  0,  6m  =  6A  =  0.  Under  these 
conditions,  solutions  can  easily  be  obtained  by  exactly  the 
same  methods  used  before,  and  they  are: 


6k 


X 

/  dx' 
o 


U(x ' ,  t  -  *12!) 


(6.1) 


and 


JT  =  -  [U(x,t)  -  U(0,  t  -  x/c  )] 

o  g  s 

+  — —  JX  dx'  L-  U (x '  ,  t  -  £—-) 

V  i  * 

- —  JX  (X-X')  dx’  —  U(x',  t  -  *1*1)  (6.2) 

4c  3  o  ^t2  cg 

o 

These  are  very  similar  to  the  solutions  obtained  before,  for 
the  initial  value  problem. 

For  the  special  case  of  U  =  U(x  -  Ct) ,  and  at  resonance, 

where  C  =  c  ,  the  solutions  reduce  to 

g 


5 
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6k 


3U 

3t 


and 


6A  =  3_  3U _ 1_  2  3^U 

K  2c  2  ^  8c  3  3t2 

g  g 


(6.3) 


(6.4) 


A  special  case  of  interest  is  evidently  that  of  a  sinus¬ 
oidal  internal  wave 


U  =  UQ  sin  K(x-Ct) 0 (t)  (6.5) 

In  this  situation,  (6.4),  for  example,  simplifies  to 

I'y 

— Zn"  UD  cos  K(x-Ct) 

Uo  sin  K(x-Ct)  (6.6) 

The  amplitude  enhancement  thus  oscillates  with  the  in¬ 
ternal  wave,  with  an  envelope  given  by 

(*t)Max  ’  *  M  ^  (6’7) 

where  x  is  the  distance  from  the  point  at  which  the  surface 
wave  is  generated.  The  growth  is  thus  linear  in  distance  for 
small  x  (i.e.  where  Kx  <<  12)  and  becomes  quadratic  in  x  for 
larger  x  (when  Kx  >>  12).  The  region  of  linear  behavior 
corresponds  to  the  "short-time"  situation  in  the  initial  value 
problem;  for  larger  Kx  the  situation  is  the  same  as  that  in 
which  the  correction  term  derived  in  Eq.  (3.7)  is  required. 
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The  result  (6.7)  is  compared,  in  Fig.  1,  with  an  ex¬ 
perimental  measurement.*  The  agreement  is  evidently  excellent, 
and  remarkable  in  that  our  perturbation  theory  result  holds  even 


VII.  VISCOUS  DAMPING 


Viscosity  is  most  easily  introduced  into  the  energy  con¬ 
servation  equation,  and  it  will  appear  as  an  additional  constant 
term  o  on  the  right-hand  side  of  Eq.  4.6,  where  1/a  is  the 
viscous  damping  time.*  Defining  a  new  amplitude  A*  by 

A  =  A' e'0t 


removes  the  term  a,  so  that  A'  satisfies  Eq.  4.6  as  written, 
without  the  o.  Hence  the  solution  for  6A'  is  given  in  3.9,  and 
the  solution  for  A,  including  viscosity,  is  simply 

‘  L  (^retarded  '  L  (Retarded]6’" 


[n  particular,  for  U  =  U(x  -  Ct)  and  at  resonance,  we  have 


6A  /  U  +  3U 

7T  =  \JT  '  1 
o  vc 


1  ce  3zl 


(7.2) 


9xz  / 


in  conformity  with  Eq.  (5.2). 


Phillips,  p.  39,  gives  numerical  values  for  o. 


r 


